Josephson Current through a Semiconductor Nanowire: effect of strong spin-orbit 

coupling and Zeeman splitting 



Meng Cheng 1 ' 2,3 and Roman M. Lutchyn 1 

1 Station Q, Microsoft Research, Santa Barbara, CA 93106-6105 
2 Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742, USA 
3 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA, 93106 

(Dated: October 26, 2012) 

We study coherent transport through a semiconductor nanowire in the presence of spin-orbit 
coupling and Zeeman splitting due to an applied magnetic field. By employing analytical and nu- 
merical techniques we develop a theory for the Josephson effect in the superconductor-semiconductor 
nanowire-superconductor structure. We show that Josephson current through the clean semicon- 
ductor nanowire exhibits a number of interesting features due to the interplay between the Zeeman 
splitting and spin-orbit coupling. We also study effect how disorder in the nanowire affects Andreev 
bound- state energy spectrum and calculate local density of states at the junction. 
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I. INTRODUCTION 

Josephson effect is an example of macroscopic quan- 
tum phenomena and is of fundamental importance in 
condensed matter physics. The applications of the 
Josephson effect are both numerous and profouncP^. 
Recently, there has been a growing interest in under- 
standing coherent transport phenomena in hybrid het- 
erostructures involving s-wave superconductors and non- 
superconducting nanostructuresSHEU ^he possibility to 
control various parameters of the nanostructures allows 
one to explore different aspects of the Josephson effect in 
these devices. 

In this paper, we study Josephson effect in 
superconductor-semiconductor nanowire-superconductor 
heterostructure (SC-SM-SC) shown in Fig [TJl. The sys- 
tem consists of a semiconductor nanowire with strong 
Rashba spin-orbit coupling(SOC) subject to a large in- 
plane magnetic field. The wire is connected to the su- 
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perconducting leads as shown in Fig. 
ical potential is in the gap, see Fig. 
one Fermi surface 12 . Coherent transport through the 
nanowire in the helical phase, which originates from 
the non-trivial interplay between large Zeeman splitting 
and strong Rashba SOC, have not been previousl y e x- 
plored^. Previous studies of the Josephson effectPHHI 
considered the case of either small Zeeman field or small 
SOC whereas we are interested in a non-perturbative 
limit where both Zeeman and SO terms are large lead- 
ing to a significant change of the nanowire bandstructure. 
Such heterostructures involving semiconductor nanowires 
and s-wave superconductors exhibit a number of inter- 
esting effects and, in particular, are very promising for 
realizing topological superconducting phase hosting Ma- 
jorana zero-energy modes^ 2 ^. We note that the present 
setup is different from the one in Ref. [20] leading to 4tt ac 
Josephson effect because the SC parts are topologically 
trivial and there are no Majorana bound states. 

In this paper, we develop a quantitative theory for 



the coherent transport through the nanowire in the 
clean limit by employing a combination of analytical 
and numerical techniques. The supercurrent through 
the nanowire in the short-constriction limit 24 25 is mainly 
transmitted by the discrete electron-hole states (Andreev 
bound states) confined in the nanowire. These bound 
states appear because an electron incident on the semi- 
conductor/superconductor interface is Andreev reflected 
as a hole with the opposite spin and momentum. To un- 
derstand how Zeeman splitting and spin-orbit coupling 
affect Andreev bound states, it is instructive to consider 
two limiting cases: (a) magnetic-field-dominated regime 
where Zeeman splitting V z is much larger than the spin- 
orbit coupling E so . In this regime, propagating electron 
and hole modes are spin-polarized, and, thus, Andreev 
reflection is suppressed, (b) spin-orbit dominated regime 
^so ^>V Z . In this case, the propagating electron and hole 
have helical character and Andreev bound states can be 
formed by the helical electrons. The intuition developed 
for the these two limits allows one to understand the gen- 
eral case of both strong Zeeman and SO couplings. We 
first consider short-junction limit and show that the he- 
lical phase in the semiconductor nanowires can be iden- 
tified by a sharp decrease of the supercurrent as one in- 
creases an in-plane magnetic field. We also consider here 
long-junction limit and show that our findings persists in 
this case as well. 

To have a better understanding of the Josephson trans- 
port in realistic situations, we consider the Andreev spec- 
trum in disordered nanowires and calculate local density 
of states (LDoS) which can be probed in tunneling ex- 
periments. It is generally found that disorder smears 
out the peaks in the local density of states correspond- 
ing to the Andreev bound states and reduces transmis- 
sion probability. However, we also show that magnetic 
field induced zero-energy crossings generically appearing 
in Andreev spectrum are robust against moderate dis- 
order. This phenomenon appears when Andreev energy 
levels are non-degenerate due to the combination of spin- 
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orbit and Zeeman couplings and is a consequence of the 
particle-hole symmetry which constraints certain matrix 
elements between a pair of zero-energy states at an iso- 
lated crossing. Thus, we predict that in this regime there 
will be generically zero-bias peaks in tunneling conduc- 
tance as a function of the superconducting phase differ- 
ence across the junction. 

The paper is organized as follows: In Sec. I we present 
the theoretical model of the superconductor-nanowire- 
superconductor junction and the scattering matrix for- 
malism for calculating the Josephson current. In Sec. II 
we examine in great detail the Andreev spectrum and 
Josephson current in short-junction limit. In Sec. Ill 
we study Josephson transport by numerically diagonal- 
izing the Hamiltonian in a finite-size system which goes 
beyond the short-junction limit. In Sec. IV we consider 
the disorder effect. 
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FIG. 1: (a) Schematic view of the SC-SM-SC junction. Semi- 
conducting nanowire is coupled to s-wave superconducting 
electrodes. The chemical potential in the nanowire can be 
tuned by the gate (not shown), (b) Energy spectrum of the 
semiconductor with strong Rashba spin-orbit coupling and 
Zeeman splitting. The two values of the chemical potential at 
fi = and \i > V z correspond to helical and non-helical states 
of the nanowire, respectively. 



II. THEORETICAL MODEL 

We consider SC-SM-SC junction shown in Fig.[T]where 
the single-channel semiconductor nanowire is modeled 
using the following Hamiltonian: 

Ho = ^(x)(- -^- f j f +V z a x -iaa y d x J^(x), (1) 

where m*,/i,a are the effective mass, chemical potential 
and strength of Rashba spin-orbit coupling, respectively; 
i/j is the spinor (^, ^|) T an d we set H = 1. The magnetic 
field B applied along the nanowire leads to the Zeeman 
energy V z = gsM^BB/2 where gsM is the ^-factor in the 
semiconductor which is assumed to be much larger than 
in the superconductor. One notices that the Hamiltonian 
is invariant under the combined operation U = icr x P of 
spatial inversion P : Pf(x) —> f(—x) and spin 7r-rotation 
ia x . The eigenvalues of Hsm in momentum space are 

s ± (k) = ^-v±Va 2 k 2 + V?. (2) 

Taking \i = to be mid-way in the gap as a reference 
point, one can see that depending on whether exceeds 
\V Z \ or not, the semiconductor can have one or two Fermi 
surfaces. The phase with only one Fermi surface is called 
"helical" because the motion of an electron and its spin 
are correlated. Conversely, if \V Z \ < /i, there are two 
Fermi surfaces and we refer to this phase as non-helical 
(i.e. normal), see Fig. |l(b)| 

Two superconducting leads on both sides of the 
nanowire are described by s-wave BCS Hamiltonian with 
Ar/l — Aoe l( P R / L being the corresponding order param- 
eters, see FigJTJa). We assume that one can neglect Zee- 
man splitting in the superconductor due to the large dif- 
ference in the Lande g- factors, e.g. the ^-factor in InSb 
is |<7sm| °° 50 whereas in the superconductor gsc ~ 2. 
Thus, one can tune the field to the regime where magnetic 



field induces large Zeeman splitting in the semiconduc- 
tor without substantially suppressing superconductivity 
in the leads. 

The quasiparticle spectrum in this system can be found 
by solving the Bogoliubov-de Gennes(BdG) equation: 

Here ^ = (u^,u±,v±,—v^) T is the wave function in the 
Nambu spinor space. 

The dc Josephson current in the junction, as a func- 
tion of the superconducting phase difference (/?, is given 
by = -^^tanh^df^e) with ^ bdng 

positive energy eigenvalues of Andreev states. Interest- 
ingly, the symmetry U defined above implies that An- 
dreev spectrum should be symmetric under cp — » — 
provided the two interfaces are mirror symmetric, which 
is assumed in our calculations. If magnetic field and the 
direction of the spin-orbit coupling are not perpendicu- 
lar to each other (i.e., there is a nonzero parallel com- 
ponent), the inversion symmetry with respect to left and 
right leads is completely broken and the spectrum can 
be asymmetric, leading to a spontaneous formation of 
supercurrent at the junction, see Refl23l for details. 

We now consider the short-junction geometry corre- 
sponding to the length of the nanowire L such that 
the dwell time of electrons r^w in the nanowire r^w ~ 
L/vp <C Aq 1 , which means that L <C £ = In 
this case, the main contribution to the Josephson cur- 
rent co mes fr om the Andreev bound states formed in the 
junctio n 1 24 ^ 25 !. The spectrum of the Andreev bound states 
in the nanowire can be obtained using scattering matrix 
formalism 2 ^. For technical reasons, we insert thin bal- 
listic normal-state leads (N) to the left and the right of 
the semiconducting nanowire allowing for a well-defined 
scattering problem 26 . The spatial separation of Andreev 
and normal scattering simplifies the problem substan- 
tially and allows one to relate the Josephson current to 
normal- state scattering matrix for N-SM-N system. 
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The spectrum of Andreev states is determined by the 
following condition^ 



det[l - r he S e (e)r eh S h (e)] = 0, 



where 



icp/2 



-i(p/2 



(4) 



(5) 



Hie 



: 7^A, r eh = 7r A ,7 : 



and 5 e (h) is the normal- state scattering matrix for elec- 
trons (holes). The particle-hole symmetry relates Sh 
to S e as Sh(s) = &S e (— e)0 _1 , where = i<j y K is 
the particle- hole operator for spin- 1/2 electrons and K 
is the operator of complex conjugation 14 . We further 
neglect the energy dependence of scattering matrices 
S e (e) ~ $ e (0), which is justified in the short-junction 
limit. We will examine the validity of the short-junction 
approximation later. We assume perfect transparency at 
the superconductor-nanowire interface unless stated oth- 
erwise. 



III. SHORT-JUNCTION RESULTS 
A. Zeeman- field dominated regime 

We consider first the case of a large Zeeman splitting 
V z ^> E so = m*a 2 , and show that the Josephson cur- 
rent exhibits two remarkable features depending on the 
value of the chemical potential: in the non-helical regime 
V z < /i, supercurrent exhibits oscillations with the ap- 
plied magnetic field whereas in the helical phase V z > \i 
the supercurrent is significantly suppressed. 

The normal electron scattering matrix can be found 
exactly in the a = limit. 



/-iB t /Al 



\ 








l/A 




1/Al 



-iB t /A* 
1/Al -iBjAl) 



hp 
sin krrL. 



sin krjL 



where the coefficients A a ,B a are defined as: 

i ( kp , k, 

k a 
kp 



cos k a L 

1 (hp 

2 L" 



(6) 



(7) 



Here kp and k a — ^2m(n — <jV z ) are the Fermi momen- 
tum in the normal leads and electron momentum in the 
nanowire, respectively. 

Having particle and hole scattering matrices, the spec- 
trum of Andreev levels can be obtained using Eq.(5) of 
the main text. In this case of V z ^ 0, a = 0, the spectrum 
is determined by the following equation: 



2 arccos 



A 



B^B\ + cos ip A^ 
arccos — — — ±arg— — +2n7r. (8) 

\ A t A i\ A i 



Here n G Z. 

Let us now consider the limit V z <C /i. In this case, the 
spectrum of Andreev levels to the leading order in V z /fi 

becomes e(ip) ~ ±Aq cos db (po^j with ipo = V z L/vp, 

and the Josephson current is given by 



2eA 



cos <po sin ^ < (p < 7r — 2cpo 
- sin (po cos ^ 7r — 2y?o < < ^ + 2<£o 
• cos cpo sin ^ 7r + 2(po < cp < 2ir 



(9) 

The current I(cp) exhibits a jump at tp = tt ± 2y?o? cor- 
responding to the zero-energy crossings in the Andreev 
spectrum. The energy spectrum and the position of the 
crossings can be understood from the doubly-degenerate 
Andreev spectrum without any Zeeman field e((p) = 
±Ao cos((/?/2), with zero-energy crossing at (p = tt. The 
Zeeman field splits the degeneracy and the zero-energy 
crossings move away from it by an amount cpo . Note that 
zero-energy crossings are not robust and appear due to 
our assumption of perfect interface transparency T = 1. 
In realistic conditions 1, and zero-energy crossings 
will become avoided crossings corresponding to the local 
minima of the energy spectrum and therefore vanishing 
of the Josephson current, which are accompanied with 
the sign reversal of the supercurrent 15 , see Figj2j 

The physical origin of the phase shift cpo can be under- 
stood quasi-classically from the Bohr-Sommerfeld quan- 
tization condition. Indeed, Andreev bound states can 
only form when the total phase acquired by the elec- 
tron undergoing two Andreev reflections is quantized: 
§ pdx + O = 27rm, m G Z. Here ^ = ±cp — 2 arccos ^ is 
the total scattering phase acquired by the electron (hole) 
at the interfaces and § pdx = =b(kf — k±)L ~ ±2V z L/vf- 
The latter agrees perfectly with the scattering matrix 
calculation. As follows from this analysis, the Joseph- 
son current in the limit of V z <C \i should oscillate with 
the Zeeman field V z with period which is proportional to 
1/L. Similar oscillations of the critic al curr ent have been 
found in SC-ferromagnet-SC junctionals. 

In the polarized limit (V z ^> /i), the propagation of 
the spin-up electrons in the nanowire is exponentially 
suppressed. The calculation of the Josephson current in 
this limit yields the following result: 

I(<p) « 2eA e- Ki (l - ^ sin 2 k x L) sin <p, (10) 



where n = y^2m*(V z — /Tj and k^ = ^/2m*(V z + ji). 
Thus, if the length of the wire £ ^> L ^> ft -1 , the super- 
current through the nanowire will be significantly sup- 
pressed. This result is strictly speaking only valid when 
a = as the presence of spin-orbit coupling leads to spin 
flips. However, the probability for the spin flips to oc- 
cur within the distance n~ l is small because E so <C V z . 
Thus, the suppression of the supercurrent in the helical 
phase should be observable even in the presence of the 
small SOC as we explicitly demonstrate below. 
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B. SOC dominated regime 



General case 



We now discuss Josephson effect in the case of small 
Zeeman splitting V z <C ^ so - The effect of magnetic field 
is perturbative here and we begin with V z = limit. 
Without any loss of generality, we can choose the direc- 
tion of the spin-orbit coupling to be along z-axis. The 
s-wave superconductor is spin-SU(2) invariant. Thus, the 
Hamiltonian for the nanowire reads 



H 



P 

2m* 



fi - apa z 



(11) 



The energy spectrum for spin- up (-down) electrons is 
shifted to the right (left). Assuming that the chemical 
potential \i > 0, we can write down the wavefunction in 
the nanowire as 



We now turn to the general case when both SO and 
Zeeman terms are large. We can develop some intuition 
by considering Andreev reflection in terms of the eigen- 
states |=b) of the Hamiltonian ([I]). From the discussion 
of the Zeeman-dominated regime, we recall that there 
are two major consequences caused by Zeeman polariza- 
tion: (a) the oscillation of the critical current with V z 
and L. (b) Significant suppression of the critical current 
when V z > fi. However, the spin-orbit interaction causes 
precession of the spin of electrons/holes within the char- 
acteristic length scale / so , which tends to wash out the 
spin polarization of electrons and holes. To this end, it 
is instructive to rewrite the s-wave pairing in the chiral 
basis: 



E c<fc<F 



b„ e 



ik(j 



)k>- (12) A^Jk^U-k) 



Here k a are defined as 
k<j = fco + rn*aa, ko 



\/2ra*/i + (m*a) 2 



Matching the boundary conditions at x — 
find the normal state scattering matrix to be 



S e = 





/ iB 





g — im* aL 


\ 


1 





iB 





gim* aL 


A* 


gim* aL 





iB 







\ o 


g—im*aL 





iB J 



(13) 
±L/2, we 



(14) 



Here 



A — cos koL 



B 



i 

~ 2 
kp 



kp 

ko 
kp 
ko 



ko 
kp 



sin koL 



(15) 



sin koL 



and kp is the Fermi momentum in the normal leads. 
The spectrum of Andreev levels follow directly: 



e((f) = ±A \ 



sin 2 ((/V2) 
1 + A 2 sin 2 k L 



(16) 



Here A = 



££) and k = ^2m*n + m* 2 a 2 , k F 
is the Fermi momentum in the superconductor. An- 
dreev levels are doubly degenerate due to the presence 
of the time-reversal symmetry as well as the approxima- 
tion S e (e) w S e (0) used here. The small corrections to 
the spectrum due to the presence of the spin-orbit cou- 
pling were considered in Refs. |18)T9] . The main effect 
of spin-orbit coupling here is actually the modification 
of the Fermi momenta which affects normal reflection 
and thus the Josephson current due to the mismatch of 
the Fermi momenta in the nanowire and superconductor. 
For special values of L satisfying the resonance condition 
koL = n7r, n G Z, transmission of electrons through the 
spin-orbit coupled wire becomes perfect and the super- 
current is the same as in the absence of SOC. 



x,y=± 



X\Ak)il>Uk)il>U-k), (17) 



where Xxx(k) = cakj \JV Z + a 2 k 2 and x~\ (&) = 

V z j \/V 2 + a 2 k 2 . Thus, both inter-band and intra-band 
pairing are generated in the chiral basis. Therefore, one 
expects that Andreev scattering at the N-S boundary has 
qualitatively similar behavior, i.e. there is both inter- 
and intra-band scattering. We can understand our ana- 
lytical results in the two limiting cases as follows: first 
in the Zeeman-dominated regime inter-band scattering 
dominates over intra-band scattering. Once the |+) band 
is gapped, quasiparticles in this band cannot propagate 
in the junction whose length is much larger than the 
length scale associated with the Zeeman gap. Therefore, 
the supercurrent will be suppressed, similar to the spin- 
polarized case discussed above. In the opposite regime 
^so ^> ^z, the Andreev reflection is dominated by the 
intra-band scattering and the situation is similar to the 
SOC-dominated regime where there is no abrupt change 
in the Josephson current across the boundary \i = V z . 

The discussion above provides some intuition about 
Andreev reflection. However, in order to understand 
the coherent transport quantitatively, one needs to con- 
sider also normal reflection at the SC-SM interface. This 
can be done by appropriately matching the boundary 
conditions and finding the corresponding scattering ma- 
trix S(e) as described in the appendix. The expres- 
sion of Josephson current in the general case is not 
particularly enlightening and we present here numeri- 
cal results instead. We proceed by first re-scaling all 
the lengths and energies in units of l so = h 2 /m*a and 
i? so = m*a 2 /h 2 , respectively. Using typical parame- 
ters for InAs m* w 0.04m e ,a w 0.1 eV • A, we estimate 
lso ~ 100 nm, Eso ~ IK. The Josephson current is mea- 
sured in units of 2eAo/h. Henceforth, all the re-scaled 
physical quantities are denoted by tilde. 

The Andreev spectra and current-phase relation for 
the SC-SM-SC junction are shown in Fig. [2j In the non- 
helical phase, Zeeman field splits the degenerate Andreev 
levels and results in zero-energy crossings. As discussed 
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FIG. 2: I(ip)/I c and E(cp) in the non-helical phase (left 
panel) and the helical phase (right panel). The parameters 
are V z = 0.1, fx = 0.5 in the non-helical phase and fx = 0.01 
in the helical phase. 



above, finite transparency of the SM-SC interface will 
lead to avoided level crossings at these points. One can 
see that as we change chemical potential from non-helical 
to helical phase, the dispersion of Andreev levels with <p 
decreases resulting in the suppression of the Josephson 
current. 

The dependence of the Josephson current on chemical 
potential and Zeeman splitting is summarized in Fig. [3] 
in which we plot \I(tt/2)\ as a measure of the size of the 
critical current. For V z ^> E so1 the suppression of the 
Josephson current between helical and non-helical phases 
manifests as a sharp borderline V z = fx in Figji^a), which 
actually persists to V z ~ E so . As one further increases 
SOC to E so ^> V z , such suppression is smeared out as 
shown in Figj3^b). The critical current oscillations with 
magnetic field are also reflected in the "stripes" in the 
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FIG. 3: Density plot of the current \I(it/2)\ as a function 
of V z and fx. In order to understand the magnitude of the 
critical current, we need to decouple the phase dependence of 
the Josephson current which is non-trivial, especially around 
<f> = 7r. We find that it is more useful to plot I {if — tt/2) 
instead of I c . In Fig. (a) we show \I(tt/2) \ for junction length 
L = 10. Fig. (b) is a zoom-in of Fig. (a) for small V z and fx 
(or equivalently, large E so ). 



region V z < fx. One can notice that as fx increases the 
onset of the oscillations requires a larger V z . Thus, the 
main features discussed qualitatively above within the 
simple models are also present in the general case. 

IV. LONG-JUNCTION RESULTS 

So far we have discussed Josephson effect for L <C £. 
We now present results for the long-junction obtained us- 
ing finite-size lattice calculation which involves full diago- 
nalization of the BdG equations. This calculation is suit- 
able for studying the Josephson effect beyond the short- 
junction limit allowing one to take into account the con- 
tribution to the current from the continuum spectrum. 
The tight-binding Hamiltonian reads 

H = ^ [fi^a'tijC^Civ' + a»ct ^(Ta^/Ci^/ + h.C.] 

+^2(V zz a-^)cl^+Y} A 4Ai + h.c). (18) 

i i 

Here we allow V^,/i, a, A to be position-dependent to 
model the junction, r are the Pauli matrices in spin 
space. We take hopping amplitudes to be Uj = 
—t(Sij-i + We assume that = a and V z \ = V z 

in the wire and zero otherwise; the SC pairing potential 
Ai is chosen to be non-zero outside of the wire with the 
magnitude Ao = O.Obt which corresponds to the coher- 
ent length £ ~ 10a. In the numerical calculations, we 
take the length of the superconductors on both sides to 
be Lsc = 200a and measure the energies and length in 
units of t and a, respectively. 

Having the excitation spectrum of the system as a func- 
tion of the phase difference one can compute Joseph- 
son current through the nanowire. The dependence of the 
Josephson current on fx and V z is shown in Fig. [4] The 
suppression of the Josephson current as one sweeps pa- 
rameters between helical and non-helical phases (fx ~ V z ) 
is also visible in the long-junction limit for small values of 
the SOC, see results for a = 0.5 in Fig. |4j As SOC is in- 
creased to a = 1.5, the sharp suppression of the current 
is smoothen out and the difference in the critical cur- 
rent between the helical and non-helical phase becomes 
less significant. Overall, we find that the interesting fea- 
tures in the Josephson effect discussed above persist in 
the long-junction limit. 

V. LOCAL DENSITY OF STATES AT THE 
JUNCTION 

In addition to measuring the Josephson current dis- 
cussed in the previous sections, one can also extract in- 
formation about Andreev energy spectrum using tunnel- 
ing spectroscopy experiments, effectively measuring local 
Density of States (LDoS). In these experiments, it is im- 
portant to take into account effect of disorder since it will 
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FIG. 4: Dependence of I c = |7(7r/2)| on \i and V z from finite- 
size lattice diagonalization. The blue(solid) lines are results 
for a = 0.5 and the red (dashed) lines are results for a = 1.5. 
The parameters for each subfigure are: (a) \i — 0.5, L = 20. 
(b) n = 0.5, L = 60. (c) V z = 0.25, L = 20. (d) V z = 0.25, L = 
60. Here the the current is normalized by the corresponding 
value at V z =0. Here p is measured from the band bottom. 



lead to a reduction of the transmission probability and, 
thus, will strongly affect the energy spectrum. For ex- 
ample, in the SNS system in the short junction limit, the 
zero-energy crossing at ip — tt as evident from Andreev 
spectrum e(cp) = ±A cos ^ is removed once we introduce 
impurities which lead to the reduction of the junction 
transparency^. The situation is different in semiconduc- 
tor nanowires in a magnetic field as we explain below. 
Let's consider LDoS defined as 



[\u na (r)\ 2 S(uj - E n ) + | (r)\ 2 S(w + E n 



(19) 

Here u a ,v a are eigenstates of the BdG equation corre- 
sponding to eigenenergy E n . In our numerical calcula- 
tions, we use Lorentzian function instead of the 5 func- 
tion and introduce some broadening of the energy levels. 
The LDoS in the clean system is shown in Fig 5^i) and 
b). As discussed above, magnetic field generically leads 
to the appearance of the zero-energy crossings in An- 
dreev spectrum. However, unlike in SNS junction, these 
crossings are robust against disorder which is a distinct 
signature of semiconductor nanowires in a magnetic field 
where Andreev states are non-degenerate in contrast with 
the conventional normal metals having spin-degenerate 
spectrum. These differences can be seen in LDoS mea- 
surements in a magnetic field which can be done exper- 
imentally since LDoS is directly related to the local dif- 
ferential tunneling conductance: 



d/(r,eV) 
dV 



oc p(r, uj — eV). 



(20) 



We now discuss the disorder- averaged LDoS p(r,u) 
where these features are visible. The disorder can be 
included in the lattice calculations by adding an on-site 



random chemical potential to the Hamiltonian ( p"8| ): 

^disorder = ^ V^cJq. (21) 

i 

The random potential Vi = V(ri) satisfies Gaussian dis- 
tribution with the correlation function (V(r)V(r , ))^i S = 
W 2 S(r — r'). In conventional metals, the strength of the 
disorder potential determines the electron mean-free path 
/ = vpT with t -1 = 2kW 2 v(ef) which is the same for 
both spins. Here v(£f) is the normal-state density of 

states at Fermi level, v(sf) = ^(^jf^)~ 1 \k=k F - I n ^ ne 
presence of spin-orbit coupling and Zeeman field, elec- 
trons in different bands have different mean-free paths 
due to, for example, Fermi velocity mismatch. In addi- 
tion, the singularity of the DoS at the bottom of the band 
leads to strong energy dependence of the mean- free path. 
Thus, unlike in conventional metals I is not a good quan- 
tity to characterize the strength of the disorder. However, 
for the lack of a better quantity to characterize the dis- 
order, we roughly estimate I w 5 for the parameters used 
in the numerical calculation and away from the bottom 
of the band. Thus, in our simulations I is much smaller 
than the length of the junction L = 20 and, thus, effect 
of the disorder is quite important here. 




FIG. 5: LDOS measured in the middle of the junction. In 
both plots, A 0.05, L = 20, a = 0.5, fi 0.5. (a) V z 
0,W = (b) V z = 0.4, W = (c) V z = 0,W = 0.3. (d) 
V z 0.4, W = 0.3. 

For the sake of concreteness, we now consider LDoS in 
the short-junction limit in the presence of spin-orbit cou- 
pling. From the scattering matrix calculation in Sec |III[ 
we already know that the spectrum of Andreev bound 
states has no zero-energy crossings, see Fig. [5|a). It is 
naturally expected that disorder does lead to any qual- 
itative changes in this case apart from broadening the 
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peaks in the LDoS and increasing the gap at ip = tt. 
This is confirmed by the numerical simulations shown 
in Fig. [5|a) and (c). The Zeeman field splits the spin 
§U(2) degeneracy and results in a complicated spectrum 
as demonstrated in FigjHJb) for the clean case and Fig. 
[5jd) for the disordered case. One can notice that there 
are two magnetic field-induced zero-energy crossings in 
the spectrum of Andreev bound states on the interval 
4> = [0, 27r]. From the disorder-averaged LDoS shown in 
Fig. [5|d), we see that the crossing at zero energy per- 
sists at least for moderate disorder strength. The reason 
of this robustness is due to particle-hole symmetry as we 
explain below. 

Let's consider an isolated crossing in the energy spec- 
trum for certain superconducting phase difference ip = 
(po. In the spirit of k • p perturbation theory, one can 
analyze effect of various perturbations around this level 
crossing. The two zero-energy states \I/o and at the 
crossing are related by particle-hole symmetry, namely 
— where 5 = <J y r y in the Nambu spinor notation 
adopted here. To understand the stability of the cross- 
ing, we consider now an arbitrary local perturbation 0{x) 
(e.g. disorder potential) and calculate its matrix element 
between zero-energy states. The perturbation 0{x) has 
to satisfy particle-hole symmetry E^O(x)E = —0 T (x). 
Using this relation, one an show that the off-diagonal 
matrix element Oio = (^i\0\^o) is zero. Indeed, one 
can show that 



O 
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(*i|O|tf ) = <*o|S + 0|tf ) 
-<^|O T £t|* ) 

dx ^(x)O T (x)E^ (x) 



(22) 



Since E T = S, one finds that SOf = —SOio = 0. 
Thus, any on-site perturbation preserving particle-hole 
symmetry does not mix \I/o and states. Another 
way of proving this fact is to show that the two zero- 
energ y l evels actually correspond to different fermion 
parity^ and as such are robust as long as there are 
no other zero-energy crossings nearby. However, generic 
perturbation produces non-zero diagonal contributions 
(^o|0|^ ) = which shift the position of the 

crossing. Thus, in the subspace spanned by the zero- 
energy eigenstates, an on-site perturbation acts as \i z . 

We now apply these results to the zero-energy cross- 
ings in the Andreev bound spectrum. As discussed in 



Sec|III[ for V z < fi there might be two such crossings 
in the spectrum located symmetrically with respect to 
ip = 7r. Each crossing has a pair of zero-energy states 
since spin degeneracy is removed by the magnetic field. 
One of the relevant perturbation near the crossing point 
is the deviation of the superconducting phase from Sip. 
According to our general argument, we can expand in 
small deviation of the phase difference Aip = ip — ipoi to 
obtain: 



Here fi z is the Pauli matrix acting in the space spanned 
by the two zero-energy eigenstates. Other perturbations, 
such as on-site disorder, also take the form 



SH' = SOooVz. 



(24) 



Combining the two contributions, we find that the spec- 
trum now reads 



E = ±\aA(p + 50 o\. 



(25) 



SHy = (aAip + bAip 2 



(23) 



Therefore, the perturbation shifts the crossing point 
— > + If the perturbation is strong enough, 

then the two crossings located at different ipi are brought 
close to each other and can annihilate each other since 
these perturbations can now mix the states correspond- 
ing to different crossings and, thus, can generically open 
a gap in the spectrum. However, such a perturbation 
needs to be larger than certain threshold which explains 
why zero-energy crossings survive in the presence of a 
moderate disorder. 



VI. CONCLUSIONS AND DISCUSSIONS 

We have developed a theory for the coherent trans- 
port through semiconductor nanowire with strong spin- 
orbit coupling and Zeeman splitting. We find a number 
of remarkable features in the Josephson current through 
the nanowire originating from the interplay between two 
competing terms: spin-orbit coupling and Zeeman split- 
ting. In particular, we show that the transition between 
helical and non-helical states of the nanowire is accompa- 
nied by a large suppression of the critical current through 
the junction when Zeeman energy is large compared to 
the spin-orbit energy scale. We support these findings 
by calculating the Josephson current both in the short- 
junction limit using scattering matrix formalism and in 
the long-junction case using the numerical diagonaliza- 
tion. The suppression of the critical current can be used 
as a diagnostic of the helical state in the nanowire. We 
also show that the critical current exhibits oscillations 
with an applied magnetic field. These oscillations are 
generic and might be relevant for recent experiments on 
topological insulators/supercoductor structures^. 

We also compute local density of states and study effect 
of disorder on the Andreev energy spectrum. Without an 
applied magnetic field, the disorder- averaged local den- 
sity of states is qualitatively similar to that of the clean 
case - the disorder mainly changes the junction trans- 
parency. In the presence of a magnetic field removing 
spin-degeneracy of the Andreev energy spectrum, we find 
zero-energy crossings which are robust against moderate 
disorder. We explain this robustness by invoking general 
arguments based on particle-hole symmetry and show an- 
alytically that any on-site perturbation only shifts the 
position of an isolated zero-energy crossing but does not 
open a gap in the spectrum. 
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Appendix A: Derivation of the scattering matrix S(s) 

We now briefly outline the derivation of the particle 
component of the scattering matrix S e (e). The system 
we consider here consists of the semiconductor nanowire 
—L/2 < x < L/2 coupled to the superconducting elec- 
trodes at x = ±L/2. The particles propagate freely 
within the nanowire —L/2 < x < L/2 in the pres- 
ence of a finite spin-orbit coupling and Zeeman field (i.e. 
a 7^ 0, V z 7^ 0). For technical reasons, we insert thin bal- 
listic normal leads to the left and right of the semicon- 
ducting nanowire^. We then calculate particle (electron) 



scattering matrix by considering N-SM-N system where 
the scattering at the SM-N interface occurs, for example, 
due to the mismatch of the Fermi momenta. The scat- 
tering at the N-SC interface is assumed to be ideal(i.e. 
normal reflection-free) consisting of Andreev reflection at 
the SC-N interface. This spatial separation allows one to 
substantially simplify the problem and relate supercur- 
rent to the normal-state scattering matrix. 

Let us first construct a basis for the particle component 
of the total scattering matrix. The wave incident on the 
nanowire corresponds to a vector: 

c in = (c+(L),c+(L),c-(R),cl(R)). (Al) 

The corresponding outgoing wave is given by 

c out = (cZ(L),c7(L),c+(R),c+(R)) (A2) 

where c^(L) correspond to right /left-moving quasiparti- 
cles with spin a in the normal region at the left bound- 
ary. The normal scattering matrix is defined as c out = 
S e (e)c in . The wavefunction in the N-SM-N system is 
given by 



' J2 4(L)^e^ x+L ^ + c-{L)^e~ ik ^ +L ^ x < -L/2 

<7=f,4. 

^(x) = I E bfXs(k s )e ik °* + b; X s(-k s )e- ik °* -L/2 < x < L/2 

S = ± 

J2 4(R)cp a e iHx - L/2) + c-(R)ip a e- ik ^ x - L ^ x > L/2 



(A3) 



Here k = ^J2m{e + /Tj where \i is the chemical potential At the right/left boundaries x = =LL/2, the wavefunc- 
in the normal metal. The spinor wavefunction ip a and tion and the probability current have to be continuous. 

For example, at x = L/2 we have the following boundary 
conditions: 



Xs(k) are given by 



cos 
sin Ok 



X+(k) 
Here tan2#fc = ak/V z . 



n = 

, X-(k) 



- sin 

cos Ok 



(A4) 



4>{L/2 + + ) =V(i/2-0+), 
d x ip(x)\ x=L/2+0 + = (d x + imaa x )tp(x)\ x= 



■ L/2-0+ 



(A5) 



Matching the boundary conditions yields a system of 
equations for c^,bf. For example, at x = —L/2 we find 



c+(L) + cf (L) = (b+e~ ik + L/2 + &;e lfc + L/2 )cos<9+ - (bU~ ik ~ L/2 - bZe ik ~ L/2 ) sin<9_ 

c+(L) + c"(L) = (b+e- ik + L/2 - b+e ik + L/2 ) sin6>+ + (bte~ ik - L / 2 + bZe ik ~ L / 2 ) cos6»_ 

k(c+(L) - cf (L)) = (k+ cos0 + + k R sm0 + )(b+e- ik + L/2 - ble ik+L/2 ) + (k R cos6>_ - fc_ sm9-)(bZe- ik - L/2 + bZe ik ~ L/2 ) 

k(c+(L) - c7(L)) = (k R cose+ + k+ sin6» + )(6+ e -' fe + L/2 + ble ik+L/2 ) + (jfe_ cos6L - k R sm9-)(bte- ik - L/2 - bZe ik ~ L/2 ) 

(A6) 
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Here fc# = ma. The equations at x = L/2 are sim- 
ilar. Eliminating bf, we obtain the particle scattering 
matrix S e (e). Using the particle-hole symmetry of the 
BdG equations one can obtain Sh(s) and compute the 
spectrum of Andreev levels using Eq.(5) of the main text. 



Appendix B: Validity of the short-junction 
approximat ion 

In the calculation of the Josephson current, we have 
made two approximations: (1) We assumed that the 
Josephson current is dominated by the Andreev bound 
states in the junction (i.e. short-junction limit), and ne- 
glected the contribution of the continuum states. (2) In 
evaluating the Andreev spectrum we have neglected the 
energy dependence of the normal state scattering matrix: 
S e (e) ~ S e (0). In this section we examine the validity of 
the second approximation. We start with the case a = 
and show that such approximation is justified apart from 
the points close to cp = 0, 2tt. 

We first derive the equation determining the Andreev 



levels with full energy dependence of the scattering ma- 
trix S e (s). When a = 0, the dispersion of the Andreev 
levels is given by 



2 arccos 
2 arccos 



e 

Ao 

€ 



arccos ^^^ ( T g)+CQ ^ + arg A ^ 



\A t (e)A±{-e)\ 

B^(-e)B±(e)+co8<p 

arccos — — — arg 

|A t (-e)^(e)| 



A t (-e)- 



(Bl) 

where A a and B a are defined in Eq. ^ with k a = 
k a (e) y/2m(ii + e- crV z ). 

Assuming that both V z and A are much smaller than /i, 
one can treat V z and e as sma ll pa rameters and perform 
perturbative expansion of Eq. (Bl ). To the leading order 



in e/fi and V z /fi, the coefficients B a and A a are given 
by: 



(JV 7 _ 



■ sin k„L 



(B2) 



Substituting into (Bl), we obtain 



2 arccos 



arccos 



(e-V z ) 2 . 



2 arccos -— = ± < arccos 
Ao I 



(e+V z ) 2 



sin 2 kpL + cos(/? 



sin kpL + cos (p 



2(s-V z )L~ 

2(e + V z )L 
v F 



} 



(B3) 



Compared to the short-junction limit, an extra term 
eL/vp appears in the equation. As a rough estimate, 
one can put in e ~ Ao. Then, this term becomes 
AqL/vf = L/£. When L/£ <C 1, we can disregard this 
term and the short-junction approximation is justified. 

However, in some cases that even for small L/£ this 
correction can result in a qualitatively different depen- 
dence on V z . Let us consider the spectrum of Andreev 
levels at ip « 0. Since the time-reversal symmetry of the 
system is broken by the Zeeman term, we expect the An- 
dreev levels to be split. However, if we start with the 
approximation of S e (e) « £ e (0), we find 



e(<p = 0) = ±A cos ^ 



V z Onk F L ± ^L 
jl v F 



which gives in the limit of V z the energy splitting 
quadratic in V z : 



Se ^ 2A ( ) kpL sin k F L. 



However, this scaling is qualitatively different if we in- 
clude the energy dependence. Indeed, assuming k F L ^> 



1, we find that 



5e, 



V z 



(B6) 



as one would expect from qualitative considerations. 
Thus, apart from the points ip « 0, 27r, the approxima- 
tion S e (e) ~ S e (0) gives qualitatively correct results. 

Let us now discuss the case with a finite SOC. Gen- 
erally, one has to rely on numerical solution but under 
the assumption that both a and V z are small with re- 
spect to /i, we obtain the following analytical equation 
for the Andreev levels, which agrees with the one derived 
in Ref. (291: 



arccos 



A 



eL 



± arcsin 



V z \sm(EL/v F ) 



E 



(B7) 

Here E = \/V z + a 2 k F . However there is a subt le p oint 



(B5) 

here: at a = the above equation reduces to rtB3|) for 



(p 7^ up to linear order in V z . However, as explained 
above Eq. (B3) around ip = contains an additional term 

(g±Vz) sin kpL 



which was neglected in Ref. [29] , 
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